arXiv:1509.03633v2 [astro-ph.SR] 25 Sep 2015 


Astronomy & Astrophysics manuscript no. ms ©ESO 2015 

September 28, 2015 


The interaction of core-coiiapse supernova ejecta with a 

companion star 

Zheng-Wei Liu\ T. M. Tauris^’^, F. K. Ropke^’^, T. J. Moriya^ M. Kruckow^ R. J. Stancliffe\ and R. G. Izzard'’^ 

' Argelander-Institut fur Astronomie, Universitat Bonn, Auf dem Hiigel 71, 53121 Bonn, Germany 
e-mail: zwliu@ynao.ac.cn 

^ Max-Planck-lnstitut fUr Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 
^ Heidelberger Institut fur Theoretische Studien, Schloss-Wolfsbrunnenweg 35, D-69118 Heidelberg, Germany 
* Zentrum fur Astronomie der Universitat Heidelberg, Institut fur Theoretische Astrophysik, Philosophenweg 12, D-69120 Heidel¬ 
berg, Germany 

^ Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 OHA, United Kingdom 
Received June 16, 2015; accepted September 10, 2015 


ABSTRACT 

Context. The progenitors of many core-collapse supernovae (CCSNe) are expected to be in binary systems. After the SN explosion 
in a binary, the companion star may suffer from mass shipping and be shock heated as a result of the impact of the SN ejecta. If the 
binary system is disrupted by the SN explosion, the companion star is ejected as a runaway star, and in some cases as a hypervelocity 
star. 

Aims. By performing a series of three-dimensional (3D) hydrodynamical simulations of the collision of SN ejecta with the companion 
star, we investigate how CCSN explosions affect their binary companion. 

Methods. We use the BEC stellar evolution code to construct the detailed companion structure at the moment of SN explosion. The 
impact of the SN blast wave on the companion star is followed by means of 3D smoothed particle hydrodynamics (SPH) simulations 
using the Stellar GADGET code. 

Results. For main-sequence (MS) companion stars, we find that the amount of removed stellar mass, the resulting impact velocity, 
and the chemical contamination of the companion that results from the impact of the SN ejecta strongly increases with decreasing 
binary separation and increasing explosion energy. Their relationship can be approximately fitted by power laws, which is consistent 
with the results obtained from impact simulations of Type la SNe. However, we find that the impact velocity is sensitive to the 
momentum profile of the outer SN ejecta and, in fact, may decrease with increasing ejecta mass, depending on the modeling of the 
ejecta. Because most companion stars to Type Ib/c CCSNe are in their MS phase at the moment of the explosion, combined with the 
strongly decaying impact effects with increasing binary separation, we argue that the majority of these SNe lead to inefficient mass 
stripping and shock heating of the companion star following the impact of the ejecta. 

Conclusions. Our simulations show that the impact effects of Type Ib/c SN ejecta on the structure of MS companion stars, and thus 
their long-term post-explosion evolution, is in general not dramatic. We find that at most 10% of their mass is lost and their resulting 
impact velocities are less than 100 kms“*. 
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1. Introduction 


The discovery of many low-mass X-ray binaries and millisecond 
pulsars in tight orbits, i.e., binary neutron stars with orbital pe¬ 
riods of less than a few hours, provides evidence for supernova 
(SN) explosions in close binaries with low-mass companions. 
The nature of the SN explosion determines whether any given 
binary system remains bound or is disrupted ( |Hills|[T98^ . An 
additional consequence of the SN explosion is that the compan¬ 
ion star is affected by the impact of the shell debris ejected from 
the exploding star ( Wheeler et al.|1975] l. Besides chemical en¬ 
richment, this kind of an impact has kinematic effects and may 
induce significant mass loss and heating of the companion star. 


Core-collapse supernovae (CCSNe) arise from massive stars. 
Determining the relationship between SNe and their progenitor 
stars remains an important problem (see |Langer|2012[ for a re¬ 
cent review). Despite the fact that the amount of SN ejecta is 
decisive for the aftermath of a SN in a binary system, the mech¬ 
anism of significant, or even complete, removal of the hydro¬ 


gen (H) envelope in H-deficient CCSNe (Type Ib/c, Ilb, and IIL) 
is not well understood. 

Although it was originally proposed that H-deficient SNe 
might arise from massive, metal-rich single stars through stellar 
wind mass loss (Woosley et al.|199^ [Maeder & Meynet||2003t 


[Eldridge & Tout||20()4| l, it has long been considered that mass 

transfer via Roche-lobe overflow (RLO) in binary systems is a 
likely mechanism to produce the progenitors of H-deficient CC- 
'odsiadlowski et al.|1992||1993[|Maund| 


SNe 


(|Paczyhski|1967[ [Pi 

et al.|2004f Eldridge et al. 2008 Stancliffe & Eldridge 2009 [ 


Yoon et al.|2blO[|Dessart et al.|20i i| l. There is growing observa- 

tional evidence that the fraction of massive stars in close binary 
systems is l arge (|Chini et al. 2012[ Duchene & Kraus 2013| l. 
[Sana et ah] ( |2012| found that more than 70% of massive stars 
are in close binary systems, which supports the idea that binary 
progenitors contribute significantly to the observed H-deficient 
SNe. 

Evidence from statistical samples of Type Ib/c SNe suggests 
that stripped, relatively low-mass, exploding stars in binaries are 
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likely to be progenitors of most Type Ib/c SNe (Drout et al. 201 ![ 


Bianco et al.| 2014] |Ly man et al. |2014[ Taddia et al.||2015| l. In 

some cases, like SN 2005ek, there is even evidence of extreme 
stripping of the SN progenitor star, with only ~ 0.1 Mq of ma¬ 
terial being ejected ( [Drout et al.|2013[ Tauris et al.|2013 2015| l. 
With HST analysis of pre-SN images, jCao et al.| ( 2013| l report the 
detection of a possible progenitor candidate at the location of a 
young, H-dehcient Type Ib SN, iPTF13bvr[^ Further light-curve 
modeling, and analysis of pre-SN and late-time photometry, of 
iPTFlSbvn suggests that an interacting binary system is likely 
to be the progenitor of iPTF13bvn rather than a standard single 
Wolf-Rayet star fBersten et al.||2014{ I Fremling et al.|[2014| El 


jdridge et al.|2015 but see Groh et al.j2013| l. All these hndings 


indicate that H-deficient CCSNe and in particular Type Ib/c SNe, 
can naturally be produced in binary progenitor systems. 

In the binary progenitor scenario of Type Ib/c CCSNe, the 
metal core of a massive star (originally 8-30 Mq) undergoes 
gravitational collapse to form a neutron star (NS), while the outer 
parts of its envelope are expelled at a high velocity and impact on 
the binary companion star. Depending on the amount of sudden 
mass loss because of the SN explosion and the momentum kick 
imparted on the newborn NS, the outcome of the SN explosion is 
either a bound system ( |Hills|198^, which may later evo lve into 
an X-ray binary system (Tauris & van den Heuvel|2006|l, or two 


single runaway stellar components, which are left behind if the 
binary is unbound ( [Tauris & Takens|1998| l. This kind of a dis¬ 
ruption of close binaries by a SN explosion has been proposed as 
a possible mechanism for producing hypervelocity stars (HVSs) 
that will leave our Galaxy because both ejected stars can receive 
high runaway velocities when asymmetric SNe are considered 
( [Tauris & Takens|1998[[Tauris|2015] l. 

After a SN explosion occurs in a binary system, the ejected 
debris is expected to expand freely for a few minutes to hours 
and eventually impact the companion star. As a result, the early 
supernova light curve can be brightened by the collision of the 


supernova ejecta with the companion star ([Kasen|2010[ 

Moriya| 

et al.|2015l[Kutsuna & Shigeyama|2015| 

Liu et al.|2015 

1 . Also, 

the companion star may he significantly 

leated and shocked by 


984[ l. Many numerical simulations of collisions of SN ejecta 
with the companion star in single-degenerate Type la SNe have 
been carried out by different groups ([Wheeler et al. 11975) [Tl 


& Fryxell[[ 19841 [Marietta et al.[[2006[ [Pakmor et al.[[20()8[ [Liu 

et al.[[2012[ [2013bt [Pan et al.[[2010[ [2()12b| l. All these studies 

suggest that about 2-30% of the companion mass, or almost the 
whole envelope of the companion star (depending on its struc¬ 
ture), is removed by the SN impact. Nondegenerate compan¬ 
ion stars are expected to survive the impact of the SN explosion 
and show some peculiar properties, such as high runaway ve¬ 
locities, overluminosity, and probably enrichment in heavy ele¬ 
ments, which may be identihed by observations. 

Whereas simulations of Type la SNe are generally performed 
based on a hxed (or small range of) explosion energy, Eej, and 
a fixed ejecta mass. Me,, CCSNe can have a large spread in Eej 


and Mq. Hirai et al. (2014i carried out two-dimensional impact 


simulations to investigate the collision of CCSN ejecta with a 
IOMq giant star companion. However, most companion stars to 
CCSNe are on the main sequence (MS) at the moment of hrst 
SN explosion in a binary system (see Kochanek|2009 and also 
our population synthesis calculations in Section[4.4[l. 


In this work, we perform impact simulations using a three- 
dimensional (3D) smoothed particle hydrodynamics (SPH) 
method to systematically study, for the hrst time, the impact of 
CCSN ejecta on MS companion stars. In Section ^ we describe 
the methods and initial models used in this work. Our impact 
simulations are presented in Section and a parameter survey 
and a detailed discussion follow in Section (4] Our conclusions 
are summarized in SectionjS] 


2. Numerical methods and models 

In this section, we introduce the codes and methods used in this 
work, and we construct the initial companion star models and 
SN ejecta models. Also, initial conditions and setup of our hy- 
drodynamical impact simulations are described in detail. 


2.1. Numerical codes 

To construct the detailed structure of the companion star at the 
moment of the CCSN explosion, we use the one-dimensional 


(ID) implicit Lagrangian stellar evolution code BEC (Danger 
1991 1998[ [Won et al.pOOb 2012[ l, which solves the hydro 


dynamic stellar structure and evolution equations ( [Kippenhahn[ 
[& Weigert[[1990|l. A recent description of the updated code is 
found in [Yoon et al.[ ( [2006[[2010| l. In our default models we as 
sume a mixing-length parameter a = l/Hp - 1.5 ( Langer|1991 
and a core convective overshooting parameter of dov = 0.10. We 
evolve stars with a given mass and solar metallicity (Z = 0.02) 
to construct ID stellar models for our subsequent impact simu¬ 
lations. 

In our hydrodynamical simulation of the collision of the SN 
ejecta on the binary companion star, we use the 3D smoothed 


particles hydrodynamics (SPH) code Stellar GADGET (Pab 
mor et al.[2012|l. An advantage of applying SPH for the problem 


the SN impact, causing the envelope of the companion star to be 
partially removed due to the stripping and ablation mechanism 
( Wheeler et al.[[T975) [Eryxell & Arnett[|l981[ [Taam & Pryxell[ 


under investigation is that no material leaves the computational 
domain, and also that momentum, energy, and angular momen¬ 
tum are strictly conserved. The GADGET code ( [Springel|2005 1 
was originally intended for cosmological simulations, hut it has 
been modihed by Pakmor et al.[(2UU8] 2012l for stellar astro¬ 
physics problems ( Pakmor et al. [20081 Liu et al.|2012[ l. As men¬ 
tioned above, several studies have successfully used the Stellar 
GADGET code to investigate the SN impact on the compan¬ 
ion star in single-degenerate SNe la (Pakmor et al.[[2008 Liu 


et al.|2012 2013c|b|al l. In particular, Pakmor et al. 120081 have 


shown that the SPH-based approach reproduces previous results 
obtained with a grid-based 2D scheme by [Marietta et al. ( 2000[ l. 
This shows that the SPH approach with the Stellar GADGET 
code captures the main dynamical effects of the collision be¬ 
tween SN ejecta and its companion star. 

2.2. Binary companion star models 

The maximum possible runaway velocities for MS HVSs ejected 
from disrupted binaries in asymmetric SN explosions have been 
investigated by [Tauris] ( [2015[ l. In their systematic Monte Carlo 
simulations of SN-induced HVSs, they chose to focus on MS 
companion stars with masses M 2 - 0.9 Mq and 3.5 Mq. This se¬ 
lection is based on observations of the main class of HVSs (G/K- 


dwarf. Brown et al. 2005 2009 2012 2014 1 and the recently dis 


covered late-type B-star HVS candidates ( Palladino et al.[: 

In this work, the same two companion stars used in Tauris[(2015|l 


* This may be the first direct detection of a progenitor of a Type Ib SN. 


are used in our impact hydrodynamical simulations. Using BEC, 
we evolved zero-age main-sequence (ZAMS) stars of these two 
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masses to an age of about 6 Myr, roughly the lower limit of their 
expected age when their massive companion star explodes. The 
structure profiles of these two companion stars at the moment of 
the SN explosion are shown in Fig. 


2.3. The supernova model 

To introduce a SN explosion, we adopt a simple analytical ex¬ 
plosion ejecta model which is constructed based on numerical 
simulations of SN explosions by Matzner & McKee ( 1999| l. We 
assume that the SN ejecta is already in homologous expansion, 
whereby the radius of a fluid element is given by r = Vej f, where 
Vej is the expansion velocity, and f is the time since the explo¬ 
sion. The density profile of the expanding SN ejecta, Pej(vej, t), 
is described by a broken power law, pg " for the outer part 
of the SN e jecta and pej or r for the inner part of the SN ejecta, 
as follows ( |Kasen|2010[|Moriya et al.|2013| l: 


Pej(Vej, 0 = < 


1 [2(5-<5)(n-5)£ej]<"-'>'" 

4n(n - 6 ) [(3 - S)(n - 3) ^ ’'"J 

1 [2(5-<S)(n-5)£ej]<'''''" _3 

4n(n - 5) [(3 - 6)(n - 3) Mej]('5-5)/2 ^ 


(Vej > V,), 

(Vej < V,), 

( 1 ) 


where Fej Mej are the SN explosion energy and SN ejecta 
mass, respectively, and 


Vt 


2(5-d)(n-5)£ej 


1/2 


(3-d)(n-3)Mej 


( 2 ) 


In SNe Ib/c, typically n ^ 10 and h ^ 1 ( [Matzner & McKee 
[1999 j l, and we adopt these to construct a homologously expand¬ 
ing SN ejecta model based on the above Eqs. (0 and (|^. Ide¬ 
ally one should start from a detailed simulation of the collapse 
of the exploding star, however, we limit ourselves to exploring 
the basic processes during the SN impact. A more realistic SN 
explosion model from hydrodynamical simulations may slightly 
change our numerical results. However, detailed parameter sur¬ 
veys are addressed in Section to mimic different SN ejecta 
models, and, therefore, we do not expect a significant change 
in our main conclusions. 

The CCSN explosion energy (Fej = 1.0 - 8.0 x 10^' erg), 
the amount of mass ejected (Mgj = 0.7 - 5.0 Mq), and the pre- 
SN orbital separation (a = 1.0 - 6.0 flmin) are treated as free 
parameters (ctmin is defined below). 


2.4. Initial setup 

Our basic setup and initial conditions for the impact simulations 
are similar to those of [Liu et~ar ( 2012 2013c|a]l. We use the 
healpix method described in Pakmor et al. ( 2012| l to map the 
profiles of density and internal energy of a ID companion star 
model to a particle distribution suitable for the SPH code (jLiuj 
jet al.|2013a| ). Before we start the actual impact simulations, the 
SPH model of each companion star is relaxed for ten dynamical 
timescales (fdyn = (Gp) *^^/2) to reduce numerical noise intro¬ 
duced by the mapping. If the relaxation succeeds, the velocities 
of the particles stay close to zero. Otherwise, we reject the SPH 
model an d redo the relaxation after adjusting the relaxation pa¬ 
rameters ( Pakmor et al.^ 2012l. We use 6 x 10^ SPH particles to 
represent the B-type star or G/K-dwarf companion star. 

The SN ejecta model described above is placed at a distance 
larger than amin from the companion star in the x-y plane. A star 



0.75 


0.69 


0.0 0.2 0.4 0.6 0.8 1.0 

M(r)/M 2 

Fig. 1. Profiles of mass density p, temperature T, and hydrogen 
abundance Ah as a function of M{r)IM 2 for a 0.9 Mq G/K-dwarf (red 
solid lines) and a 3.5 Mq late-type B-star (blue dashed lines) companion 
star model at the moment of the CCSN explosion. Here, M(r) is the 
mass enclosed within radius r, M 2 is the total mass. 


in a binary cannot be closer to its companion than the Roche- 
lobe radius because mass transfer would result. We therefore use 
the Roche-lobe radius of MS companion star ( |Eggleton|1983|l to 
set the minimum possible separation amin in our systems (jTaurisj 
[2^ , 


^min 


0.6g2/3 +^1/3) ^ 

0.49^2/3 ^2 


(3) 


where q is the mass ratio between the companion star (M 2 ) 
and the exploding star (Mhc), i.e. q - M 2 /MHe, and R 2 is 
the radius of the companion star ( jEggletori 1983| l. As a con¬ 
sequence, the minimum orbital separations between the explod¬ 
ing star and our companion stars are 2.74 Rq (G/K-dwarf com¬ 
panion) and 5.56 Rq (late-type B-star companion), correspond¬ 
ing to Amin = 3.55 R 2 and = 2.55 R 2 , respectively. We 
assume that the exploding star is an evolved helium star of a 
mass Mhc = 3.0 Mq which corresponds to a SN ejecta mass of 
ss 1.4 Mq with the formation of a NS with a gravitational mass of 
1.4 Mq, given that about 0.2 Mq is released by neutrinos from 
lost gravitational binding energy during the explosion. The im¬ 
pact of the SN ejecta on the binary companion is simulated for 
more than 6000 s, at which point the amount of mass removed 
from the companion envelope and its resulting impact velocity 
have converged (Eig. [^. 


3. Hydrodynamical results 

In this section, our impact simulations of a G/K-dwarf compan¬ 
ion star model (M 2 = 0.9 Mq) with an initial binary separation 
of a = 5.48Ro (~ 7.1 R 2 ~ 2.0 Amin) are chosen as an exam¬ 
ple. The detailed evolution of the companion star during the 


Article number, page 3 of 


12 



































































Table 1. Resolution test for a G/K-dwarf companion model. 


Nstar 

^total 

AM2 

[Mq] 

Vim 

[kms-'l 

A/Wacc 

[Mo] 

5 000 

12 640 

0.0327 

60.4 

3.7 X 10-4 

10 000 

25 427 

0.0256 

55.8 

8.1 X 10-"^ 

50 000 

127 656 

0.0123 

37.8 

8.8 X 10-"^ 

100 000 

255 513 

0.0105 

33.4 

9.3 X 10-4 

500 000 

1 278 568 

0.0066 

27.4 

9.5 X 10-4 

1 000 000 

2 556 994 

0.0062 

24.8 

9.7 X 10-4 

2 000 000 

5 113 998 

0.0057 

22.6 

9.6 X 10-4 

4 000 000 

10 227 949 

0.0055 

22.5 

9.8 X 10-4 

6 000 000 

15 341 986 

0.0055 

22.6 

9.7 X 10-4 


Notes. Nstar and Ntota\ are the number of particles used to represent the 
companion star and the binary system, respectively. All particles have 
the same mass Msph. 



Fig. 2. Amount of mass removed from the companion star as a 
function of time in the impact simulation of a G/K-dwarf companion 
(M 2 - 0.9 Mq, a = 5.48 Rq ~ 7.1 R 2 } with different resolutions, i.e., 
different total numbers of SPH particles in the simulation (Table[^. The 
standard resolution used in our simulation is highlighted in red. 


impact of the SN ejecta is qualitatively described. Under the as¬ 
sumption of homologous expansion, a typical SN ejecta model 
in our simulation is constructed according to Eqs. Q and 
with Eej = 1.0 X 10^' erg and Mgj = 1.4 Mq . The maximum 
velocity at the ejecta front is set large enough (^ 30 000km s“') 
that the assumed is actually contained in the ejecta. 


3.1. Resolution test 

To assess the reliability of our numerical results, we have per¬ 
formed a convergence test to check the sensitivity of our results 
to different resolutions by varying the particle numbers from 
5 000 to 6 X 10^ in the companion star (Table [^. Because all 
our SPH particles have the same mass, the number of SPH parti¬ 
cles representing the SN ejecta is fixed once the SPH companion 
model is created. 

The amount of removed companion mass as a function of 
time in our impact simulations on a G/K-dwarf companion star 
model (a = 5.48 Rq) with different resolutions is shown in 
Fig. 0 To calculate the amount of mass removed (evaporated) 
from the companion star by the SN impact, we sum over the 
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Fig. 3. Amount of mass removed from the companion star (top panel) 
and the impact velocity caused by the SN ejecta (bottom panel) as a 
function of time after the SN explosion in our impact simulations of a 
G/K-dwarf companion (M 2 = 0.9 Mq) with an initial binary separation 
of a = 5.48 Rq ~ 7.1 R 2 ~ 2.0 Here, a„,in is the minimum possible 
separation of the binary system. 


masses of all unbound SPH particles (stripped-i-ablated parti¬ 
cles) that originally belonged to it. Whether or not a particle 
is bound to the star is determined by calculating its total energy, 
Ftot = Fkin + Fpot + Fin, where E^m, Epot and Ein are the kinetic 
energy (positive), potential energy (negative) and internal energy 
(positive), respectively. If Ftot > 0, the particle is unbound. The 
center-of-mass motion of the star is subtracted when calculat¬ 
ing the kinetic energy of each particle. To make the difference 
between the results obtained with particle numbers above 2 mil¬ 
lion better visible, only the results with a resolution ranging from 
about 1.3 X 10^ to 1.5 X 10^ SPH particles are plotted in Fig.|^ 
The difference in removed companion mass between simulations 
with 5.1 X 10^ (the solid line) and 1.0 x 10^ (the dash-dotted 
line) SPH particles is smaller than 4%, the difference between 
1.0 X 10^ and 1.5 x 10^ (the dotted line) SPH particles is less 
than 2% (Table [T]i. We find that the amount of removed com¬ 
panion mass (AM 2 ), the impact velocity (vim) and the amount of 
accreted ejecta mass (A/Wacc) are sufficiently converged (Table[^ 
when using more than 4 x 10^ SPH particles to represent the 
companion star and we thus settle on the number 6x10^ (which 
corresponds to a total of about 10^ particles in the simulation). 
The mass of a single particle, /Wsph, is then 10“® - 10“^ Mq. 

3.2. Description of evolution during/after SN ejecta impact 

Figure [^demonstrates the typical evolution of the gas density in 
the orbital plane in our impact simulation of a 0.9 Mq G/K-dwarf 
companion model with an initial binary separation of 5.48 Rq ^ 
7.1 R 2 . The basic impact processes are quite similar to previous 
hydrodynamical simulations with a MS donor model in normal 
SNe la dLiu et al.|2012l|Pan et al.|2m2bl ). 

At the beginning of the simulation, the MS companion star 
is in equilibrium (Fig.j^). The SN explodes on the right side of 
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Fig. 4. Density distributions of all gas material as a function of the explosion time in the impact simulations for a G/K-dwarf companion model 
(M 2 = 0.9 Mo) with a binary separation of 5.48 f?o (= 7.1 R 2 )- The direction of motion of the incoming SN shell front is from right to left (see 
arrow symbols). The impact is initiated at f 100 s after the explosion. The color scale shows the logarithm of the mass density in g cm“^ . 


the companion star, the SN ejecta expands freely for a while and 
hits the MS companion star, removing H-rich material from its 
surface through stripping (momentum transfer) and ablation (SN 
shock heating). The latter is caused by a bow shock propagating 
through the star, leading to internal heating and additional loss 
of mass from the far side of the companion star (Fig. [^). In 
the final snapshot (f = 5000 s), all removed companion material 
caused by the SN impact is largely embedded in a low-velocity 
SN debris behind the companion star. At the end of the sim¬ 
ulation, we calculate the amount of removed companion mass, 
AM 2 ~ 0.0055 Mq. Mass loss caused by the SN impact reaches 
a steady state by the end of the simulation, we therefore do not 
expect any further loss (Fig.[^. In Fig.[^ a small peak of the re¬ 
moved mass is seen in the early phase after the SN ejecta hits the 
companion star. This is because a shocked particle in the center 
of the star can have a larger kinetic than potential energy when 
the shock wave passes through the star. However, it is not able to 
leave the star and soon loses its kinetic energy to some neighbor 
particles and therefore becomes bound again. 

Despite the fact that the kinetic energy of the SN ejecta inci¬ 
dent upon the companion star is much greater than the total bind¬ 
ing energy of this star, the star survives the impact because en¬ 
ergy is deposited in only a small fraction of the stellar mass and 
excess energy is redirected into kinetic energy of the expelled 
material ( |Fryxell & Arnett|1981[ |Taam & Fryxell|1984| ). How¬ 
ever, the companion star receives linear momentum and hence 
a resulting impact velocity because of the collision with the SN 
ejecta, causing the star to move in the same direction as the SN 
ejecta with a speed of 22kms“'. Furthermore, the surviving 
companion star is forced out of thermal equilibrium by internal 



Fig. 5. Mass fraction of companion material in the SN ejecta in 
our simulation of the G/K-dwarf companion star model with a binary 
separation of 5.48 7?© (= 7.1 R 2 ). The purple end of the color scale 
corresponds to pure SN ejecta material, while the red end of the color 
scale represents pure companion material. 


heating from the passing shock wave. Hence it is expected that it 
remains bloated on a thermal timescale (up to tens of Myr) after 
the SN explosion (cf. Section |4.5| l. To illustrate how the su- 
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Fig. 6. Total removed companion mass (left panel), resulting impact velocity of the companion star (middle panel) and the amount of accreted 
contamination from the SN ejecta (right panel), as a function of initial binary separations for a G/K-dwarf (square symbols, M 2 = 0.9 Mq and 
R 2 ~ 0.77 Rq) and a late-type B-star (filled circle symbols, M 2 = 3.5 Mq and R 2 ~ 2.18 Rq) companion model. The results of our example model 
(Section]^ are shown with red filled square symbols. Power-law fits are also presented in each panel. 


pernova ejecta are mixed with the removed companion material 
after the SN explosion, the amount of material that originally 
belonged to the companion star is shown relative to the total 
amount of material for the simulation of a G/K-dwarf companion 
star model in Fig. It is evident that most of the removed com¬ 
panion material is confined to the downstream region behind the 
companion star, creating a hole in the SN ejecta with an opening 
angle of ^ 30° with respect to the .r-axis in our simulation. 

Qualitatively, the impact proceeds in the same way in all our 
simulations. The numbers given above are derived in the special 
case of a simulation using Eg = 1.0 x 10^* erg. However, as 
we demonstrate systematically below, increasing the explosion 
energy enhances the effects of the SN ejecta impact. 


3.3. Impact velocity of the companion star 


The impact velocity kick (up to about 100 km s *, cf. Section]^, 
which the companion star receives as a consequence of absorp¬ 
tion of linear momentum from the SN ejecta, is in the orbital 
plane and opposite to the direction to the exploding star. Al¬ 
though this impact velocity is small compared to the pre-SN or¬ 
bital velocity of the companion star, it increases the probability 
of disrupting the binary system and, in that case, also contributes 
(although with a relatively small amount) to the escape velocity 
of the ejected star ( Tauris & Takens|1998 Tauris|20T5) l. In ad¬ 
dition, this impact velocity affects the eccentricity of the binary 
system if it survives the SN explosion. 


The symmetry-breaking effect of orbital motion is ignored 
in our impact simulations because the orbital and spin velocities 
of the companion stars are much lower than the typical ejecta 
expansion velocity of a few lO^kms^'. We expect that some 
small morphological differences, such as the post-impact density 
distribution of all removed companion material and SN ejecta, 
would be seen if the orbital motion of the companion star is in¬ 
cluded, but it should not affect the basic results presented in this 
work (see also|Liu et al.|2013b|. 
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3.4. Contamination from the SN ejecta 


After the interaction with the SN ejecta, the envelope of the 
companion star may be enriched by heavy elements from the 
SN ejecta while part of its H-rich material is stripped off by 
the SN impact. As a consequence, a surviving companion star 
may show unusual chemical signatures compared to a normal 
MS star, in case a significant amount of SN-ejecta material is ac¬ 
cumulated onto it. These chemical peculiarities, if detected, can 
potentially differentiate between SN-induced HVSs from HVSs 
ejected from the Galactic center ( Hills|19^ or from d ense stel¬ 
lar clusters by dynamical interactions ( Leonard|1991| l. Possibly 
even the explosion mechanism (thermonuclear vs CCSN) could 
be determined for these chemically enriched HVSs ( |Schaffen~ 
|roth et al.|2015| ). 


In our hydrodynamical simulations, we trace at every 
timestep all bound particles that originally belonged to the SN 
ejecta after the explosion. It takes some time (about 4000 s 
in our typical G/K-dwarf companion model) for this captured 
ejecta material to settle onto the surface of the companion star. 
At the end of the simulations, the amount of ejecta mass accreted 
onto the surface of the surviving companion star remains at a 
(roughly) constant value. In our example model (Fig.|^, we find 
a total ejecta mass of ~ 10~^ Mq is accumulated onto the surface 
of the companion star at the end of the simulation compared to 
the removal of 5.5 x 10“^ Mq of H-rich material. Unfortunately, 
we cannot provide detailed information of the chemical elements 
of the SN ejecta material that accumulate onto the surface of 
the star. This is because in our simulations only a simple ejecta 
model is constructed according to Eqs. ([^ and (|^, without any 
information about the chemical composition the ejecta material. 


From our simulations, we find that most of the contamina¬ 
tion is attributed to particles with low expansion velocity in the 


SN ejecta, i.e. the innermost region of the SN ejecta model (Liu 
|et al.||^13a| l. This is explained by the lower kinetic energy of 
those particles which makes it easier to accrete onto the surface 
of the companion star after the momentum transfer. Therefore, 
we argue that the composition of the ejecta material that pollutes 
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the companion star is sensitive to the nuclear burning at the cen¬ 


ter of the explosion. Based on CCSN modeling (e.g., Nomoto 


et al. 20061, we expect a potentially observable chemical enrich¬ 
ment of Fe, Si, and possibly Ca (e.g., Schaffenroth et al.]|2015 


Sabach & Soker 2015|l, if mixing processes (e.g., Stancliffe et al. 


2007[ STancliffe & Giebbeek||2008| l are not efficient on a long 


timescale. Simulations with a more detailed model of the ex¬ 
ploding star are needed to address the issue of chemical contam¬ 
ination in a future study. 

In addition, we point out that the SPH method suppresses hy- 
drodynamical instabilities seen in grid-based simulations, which 
may affect the estimated contamination from SN ejecta. How¬ 
ever, earlier impact simulations with the same SPH approach 
(|Liu et al.||2013c] l and 3D grid-based simulations ( |Pan et aT] 
|2012b| l f^nd the same order of magnitude in the amount of 
nickel contamination on a He star donor from SN la ejecta. 
This indicates that the degree of instability captured in our SPH 
method is sufficient for our study, or these instabilities may not 
play an important role in determining the amount of accreted SN 
ejecta in the simulation. This, however, should be tested in more 
detailed future studies. 


4. Discussion 

The kinetic energy, mass of the SN ejecta, the separation be¬ 
tween the exploding star and its companion, and the structure 
of the companion star at the time of the explosion are expected 
to be four major physical parameters of the progenitor system 
that influence the dynamics of the SN ejecta impact on the com¬ 
panion star. Different binary systems evolve to various evolu¬ 
tionary stages and have distinct binary parameters when the SN 
explodes. Consequently, some properties of the companion star 
and the binary system differ significantly from our example mod¬ 
els. For instance, the binary separation (a), the mass of the ex¬ 
ploding star (Mhc, which leads to different ejecta masses), and 
the structure (mass, M 2 , and radius, R 2 , and thus density profile) 
of the companion star at the moment of the SN explosion, are all 
expected to be different in various binary systems. 

In this section, we treat the binary separation (a) at the time 
of the explosion, the ejecta mass (Mg), or the SN explosion en¬ 
ergy (Egj), as independent physical parameters, and we carry out 
a parameter survey to investigate the dependence of our numeri¬ 
cal results on these parameters. Predictions from our population 
synthesis calculation are also discussed. 

4.1. Binary separation dependency 

To test the influence of the ratio of binary separation to the com¬ 
panion star radius, a/R 2 , we adjust the binary separation for 
a given companion star model. A fixed SN explosion model 
with an ejecta mass Mgj = 1.4 Mq and an explosion energy 
Fej = 1.0 X 10^' erg is used to model the SN event. 

Figure shows the effects of varying the orbital separation 
parameter a/R 2 , by a factor of about 6, on the total amount of 
removed companion mass (AM2), the resulting impact velocity 
(vim), and total accumulated ejecta mass (AMacc), for the 0.9 Mq 
and 3.5 Mq companion star models. It shows that all three quan¬ 
tities decrease when increasing the orbital separation of the pre- 
SN binary. For instance, AM 2 decreases by a factor of almost 10 
as the parameter a/R 2 increases by a factor of 2. The relations 
can be approximately fitted with power-law functions (Fig. |^, 
which are consistent with previous impact simulations of ther¬ 
monuclear (la) SNe (e.g., [Marietta et al.|[2000t jPakmor et al.j 


0.12 

0.10 

^ 0.08 
< 

K 0.06 

ra 

E 

"S 0.04 
o 
E 
cu 

CC 0.02 


0.0 2.0 4.0 6.0 8.0 10.0 

Explosion energy [10®^ erg] 

Fig. 7. Total removed companion mass (AM 2 , in a unit of M 2 ) as a 
function of the SN explosion energy for our 0.9 Mq G/K-dwarf com¬ 
panion model. All parameters except the SN explosion energy are kept 
constant, meaning that we only change the SN explosion energy ac¬ 
cording to Eqs. Q and Thus, for a fixed initial binary separation 
of 5.48 (2amin) and a total ejecta mass of 1.4 Mq, different SN ex¬ 
plosion energies of (1.0, 4.0, 6.0, 8.0) x 10^' erg are investigated. The 
fitting function and fitting parameters are also shown. 
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Fig. 8. As Fig.|^ but for the resulting impact velocity (vjn,) and bound 
SN ejecta mass (AMacc). 


[20081 [Liu et al.[[2012[ [Pan et al.|2012b| l and CCSNe with a red 
giant companion star ( Hirai et al.|2014 i. 

A more careful inspection of Fig.j^shows that the decrease in 
AM 2 weakens systematically as the binary separation increases 
to be large, so systems with a/R 2 ^ 8 deviate slightly from the 
power law. This indicates that the detailed interactions are possi¬ 
bly more complex. For example, the amount of evaporated com¬ 
panion star material, and the accreted (bound) SN ejecta mass, 
may be poorly determined at long orbital separation because nu- 


Article number, page 7 of 


12 







































Fig. 9. As Figs. 1^ and 1^ but as a function of the SN ejecta mass. 
All parameters (binary separation of 5.48 Rq and a SN explosion energy 
of 1.0 X 10^' erg) except the SN ejecta mass (Mej) are kept constant. 
Different values of M^j = 0.7, 1.4, 2.5, 4.0, and 5.0 Mq are applied. 



Fig. 10. As Fig.|^ but for the amount of mass accreted from the SN 
ejecta (AMacc) as a function of the SN ejecta mass (Mej). 


merical difficulties were encountered when the removed mass (or 
accreted SN ejecta mass) is very small, i.e., when only a small 
number of SPH particles are involved. Furthermore, it cannot 
be ruled out that the impact simulation process also depends on 
the time evolution of the density of the SN ejecta hitting the star. 
However, we do not expect these uncertainties to strongly affect 
our basic results and conclusions, and the power-law fits pre¬ 
sented here are fairly good approximations for the full range of 
alR2. 

4.2. SN explosion energy dependency 

The influence of the SN explosion energy on the ejecta interac¬ 
tion with the companion is studied on the basis of our G/K-dwarf 
companion star model with a ^ 5.48f?© ^ 7.1 R 2 . All parame¬ 
ters but the SN energy are kept constant. With a fixed ejecta mass 
of Mej = 1.4 Mq, we create four different SN ejecta models, us- 
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Fig. 11. Top: momentum within each velocity shell of initial SN ejecta 
as a function of ejecta velocity. Bottom: total momentum of outer SN 
ejecta within regions with ejecta velocities > 5000, 6000, 7000, and 
7500km s“'. Different ejecta masses of Mej = 0.7, 1.4, 2.5, 4.0 and 
5.0 Mq are applied. 


ing Eqs. (fTI and (ItJi, with ejection energies (1.0-8.0)xl0^’ erg 
dTaddia et al.|2015| ). 

Numerical results for the total removed companion mass, 
AM 2 , the resulting impact velocity, vim, and total accreted mass 
from the SN ejecta (i.e., contamination), AMacc, are illustrated 
in Figs. and 1^ as a function of the SN explosion energy, 

The correlationDetween AM 2 and Eej is, to a good approxima¬ 
tion, linear (Fig.|^. As expected, AM 2 increases with E^j. The 
behavior of both Vim and AMacc are better described by power- 
law functions (Fig.j^. Whereas Vim increases with Eej, AMacc is 
smaller at larger E^j. 

4.3. SN ejecta mass dependency 

To investigate the influence of the SN ejecta mass, Mg, on the 
outcome of the SN impact, (i.e., AM 2 , Vjm and AMacc) we vary 
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the SN ejecta mass by fixing all other parameters, i.e., - l.Ox 

10^^ erg, a = 5.48/?o, and apply a maximum ejecta velocity ^ 
30000kms *. Under these assumptions, five SN ejecta models 
with ejecta masses 0.7, 1.4,2.5,4.0, and 5.0 Mq are created using 
Eqs. Q and (|^. 

Figs.[9]and [T0|show the effect of the varying SN ejecta mass 
on the interaction between the SN ejecta and the companion star 
for our 0.9 Mq G/K-dwarf companion model. With our method 
for creating different SN explosion profiles, we find that the SN 
ejecta does not significantly affect AM 2 (Fig. [^). However, it 
is found that the impact velocity of the companion star (Fig. [^) 
decreases nonmonotonically with increasing ejecta mass. This is 
a somewhat surprising and perhaps counterintuitive result given 
that for a fixed total kinetic energy, Fg = 1.0 x 10^' erg, the in¬ 
cident momentum (i.e., the momentum of the SN ejecta that hits 
the star) is expected to increase by a factor of 2.7 as the ejecta 
mass increases from 0.7 to 5.0 Mq (using Fej = P^/(2Mej)). 
Therefore, the impact velocity of the star is expected to increase 
by a factor of about 2.7 by the analytic analysis of [Wheeler et al.| 
( |1975| l, which is, however, not seen in our simulations. There 
are several reasons for this discrepancy between expected val¬ 
ues from the analytic work and the results of our simulations. 
In Fig. it is shown that the impact velocity of the star almost 
reaches a maxium value about 600-700 s after the SN explo¬ 
sion. This suggests that the impact velocity of the star is mostly 
determined by the momentum transfer during the collisions with 
high-velocity ejecta at early-time phase. Late-time interactions 
between the SN ejecta and the star do not significantly increase 
the impact velocity further because the drag on the companion 
star is significantly reduced at the late-time phase ( [Fryxell & Ar-| 
|nett|198l1 l, and the momentum transfer due to ablation (which is 
dominant at the late-time phase) is quite inefficient. In addition, 
the shape of the companion star is changing while its mass is 
being removed by the SN impact, which may reduce its cross- 
sectional area and thus reduce the momentum transfer and drag 
efficiency ( [Fryxell & Arnett||I98I| l. All these effects were not 
included in the analytic calculations of [Wheeler et al.| ( |1975| l. 
Therefore, the impact velocity of the companion star is very sen¬ 
sitive to the total momentum of the very outer SN ejecta layers 
with high expansion velocities. We refer to [Fryxell & Arnett| 
( 1981| l for further detailed descriptions of the impact velocity of 
the companion star after the SN impact. 

In Figure [TT^, we present the momentum within each spher¬ 
ical velocity shell of the initial SN ejecta as a function of ejecta 
velocity for different ejecta mass models. In addition, the total 
momenta of the outer SN ejecta regions with an expansion ve¬ 
locity in excess of 5 000,6 000,7 000, and 7 500km s“' are cal¬ 
culated for different ejecta mass models, respectively (Fig.fTTjr). 
Our ejecta model leads to most of the momentum being dis¬ 
tributed in the inner part, i.e., low expansion velocity region, of 
the SN ejecta rather than in the outer part, i.e., high expansion 
velocity region, as we increase the SN ejecta mass (Fig. [TT^). 
For example, as shown in Fig. [TT| r, the total momentum within 
a region with expansion velocities larger than 7000 km s ' in¬ 
creases first, then decreases with the increase of the SN ejecta 
mass. This high-velocity part of the SN ejecta hits the star within 
< 600 s after the explosion and is expected to be the main factor 
that determines the final impact velocity of the companion star 
(Fig. s- Therefore, it is not surprising that the dependence of the 
impact velocity (vim) on the SN ejecta mass (Mej) in our simula¬ 
tions follows a pattern as shown in Fig[^, keeping in mind that 
the impact velocity of the star strongly depends on the momen¬ 
tum distribution of the outer part of the SN ejecta. Adopting a 
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Fig. 12. Population synthesis distribution of the companion star mass 
(M 2 ) as a function of the binary separation (a) in case the SN explodes 
as a Type Ib/c. Nothing is plotted in the regions with number fraction 
smaller than 10“^^. 


SN ejecta model that is different from ours is therefore expected 
to yield a change in the resulting impact velocity of the star. 

As expected, AMacc increases with larger Mg and is approxi¬ 
mately fitted by a power-law function (Fig.[T0|. As already men¬ 
tioned in Section 3.4 most of the ejected material captured by 
the companion star comes from the innermost low-velocity part 
of the SN ejecta ( Liu et aL]|2013c| l. The enclosed mass of the 
low-velocity ejecta increases as Mgj increases, which leads to 
more ejecta material being captured by the companion star after 
the momentum transfer. 


4.4. Population synthesis predictions 

To estimate the fraction of CCSNe in which light-curve bright¬ 
ening by the collision of SN ejecta with a stellar companion can 
be observed, Moriya et al. ( 2015| l perform populations synthe¬ 
sis calculations for SNe with the binary_c/nucsyn code ( jlzzardj 
jet al.|[2004j |2006| |2009| l. Under an assumption of a Galactic 
star formation rate of 0.68 - 1.45 Moyr^' and the average stel¬ 
lar mass of the Kroupa initial mass function (0.83 Mq), the total 
CCSN rate of the Galaxy predicted from their population syn¬ 
thesis is 0.93-1.99 X 10-2 yr-' (36% SNe Ib/c, 10% SNe Ilb, 
54% SNe II), consistent with the estimated Gal actic CCSN rat e 
of 2.30 + 0.48 X 10-2 yr“* from recent surveys ( jUi et al.|2011 1 . 
In this section, we use the data from their population synthesis 
calculations directly to discuss the effects of an explosion on the 
companion star in SNe Ib/c. The basic parameters and assump¬ 
tions in their calculations are summarized as follows. The initial 
mass function of [Kroup^ ( |200 1 j ) is applied in combination with a 
flat distribution of mass ratios (q - M 2 /Mi) and the assumption 
of circular binary orbits. The stars have a metallicity Z = 0.02. 
The orbital period distribution follows that of [Sana et al.|(|2012|l 


for O-type primary stars in excess of 16 Mq, Raghavan et al. 
pOlOj l for solar-like primary stars with a mass below 1.15 Mq 
and a fitting function interpolated linearly in primary mass be¬ 
tween these two limiting masses. The common-envelope ejec¬ 
tion efficiency is set to qqe = 1 and the envelope binding energy 
parameter. A, is adapted from ( jPewi & Tauris|2000| [200 1| [Tau-j 
ris & Dewi|2001|l. Systems with primary masses 3 - 80 Mq and 


mass ratios between 0.1 Mq/Mi and 1.0 are evolved with orbital 
periods between 0.1 and lO*" days. 
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Fig. 13. Density distribution simulations of the companion star using a 3.5 Mq B-star MS model with different initial binary separations. Only 
the bound material that originally belonged to the companion star and final accreted ejecta are shown. The logarithm of density is color-coded. 
The top-left panel shows the initial configuration for all three cases. The motion of the incoming SN ejecta is from right to left. The following 
three panels (for initial orbital separations of = 2.55, 5.10 and 10.20, respectively) show snapshots of the mass density at the moment that the 
removed companion mass, the resulting impact velocity, and the contamination from the SN ejecta, have converged to (almost) constant values. 
In the widest system (lower-left panel), the impacted star is not significantly inflated compared to its initial state before the SN impact. The white 
curves show constant density contours (logiop = -2.0). 


Figure [represents the distributions of the companion mass 
and the ratio of binary separation to companion radius (a/f? 2 ) at 
the moment of the explosion for SNe Ib/c. Most SNe Ib/c have 
an orbital separation of > 5.0 R 2 , which is about a fraction of 
> 95% in our binary population synthesis calculations. Further¬ 
more, with the distributions of a //?2 in Fig- 12 we can simply 
estimate AM 2 , Vim and AMacc by applying our power-law rela¬ 
tionships stated in Fig.j^ Typically, we find that less than 5% of 
the companion mass is removed by the Type Ib/c SN impact and 
the star receives an impact velocity of only a few 10km s *. In 
addition, typically a few 10“^ Mq of SN ejecta mass is found to 
accumulate onto its stellar surface after the SN impact. 

The above results are rough estimates based on the power- 
law relationships obtained from our simulations for a G/K-dwarf 
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and a late-type B companion star model (Fig. j^. However, 
different fitting parameters are found for these two companion 
stars, which indicates that the detailed structures of the compan¬ 
ion star could affect our numerical results. Based on previous 
simulations ( |Liu et al.|2012| l the fitting parameters are expected 
to change slightly for different companion structures. Neverthe¬ 
less, we do not expect significant discrepancies because, accord¬ 
ing to our population synthesis, a major fraction of the compan¬ 
ion stars (> 85%) are found to be in their MS phase when the 
Type Ib/c SNe occur. 
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4.5. Surviving companion stars 


After the SN explosion, the binary system may be disrupted. As 
mentioned in Sections |3.2| and |3.4| the companion stars surviv¬ 
ing the SN explosion probably exhibit some peculiar luminosity 
and chemical signatures and may thus be identihed (e.g., |Pan| 
|et al.|2012al l if these companion stars are strongly impacted and 
significantly heated by the SN ejecta. 

Figure [1^ shows the mass density structures of surviving 
companion stars from our simulations of a late-type B-star com¬ 
panion model with different initial orbital separations. Again, 
we use the default SN ejecta model described in Section |4.1| 
At small orbital separations, the companion star is strongly im¬ 
pacted and heated by the passing shock wave and hence it is 
bloated after the SN impact (top-right panel of Fig. [T3| l. As the 
orbital separation increases to ~5 R 2 , the SN impact and heating 
of the companion star are reduced. Only a small amount of com¬ 
panion mass is removed and the SN heating is too inefficient to 
make the surviving companion star inflate by much (bottom row 
of Fig. [^. 

Our population synthesis modeling suggests that most CC- 
SNe (> 95%) occur in systems with an orbital separation larger 
than about 5 .QR 2 (Section [4.4| i. We caution that there are large 
uncertainties in population synthesis studies, which may influ¬ 
ence the results. In this case, these mainly relate to the in¬ 
put physics of common-envelope evolution and the subsequent 
Case BB RLO from the naked helium star prior to its explo¬ 
sion ( |Tauris & van den Heuvel||2006| and references therein). 
In addition. Fig. 12 shows that at the moment of SN explosion 
most companion stars have a higher mass than our 3.5 Mq late- 
type B-star companion. More massive stars have higher surface 
escape velocities and hence less companion mass is removed. 
Also, the fractional change in the internal stellar energy created 
by the shock wave is smaller because of higher internal energy 
in more massive stars. Consequently, if our population synthesis 
is correct, we do not expect that SN explosions significantly af¬ 
fect the detailed structures of companion stars in most CCSNe. 
This also implies that the effect on the post-impact evolution of 
the companion stars is not expected to be very dramatic. There¬ 
fore, some observable signatures, e.g., swelling and overlumi¬ 
nosity, predicted for MS companions in SNe la systems are, in 
general, not expected to be seen dramatically to the same extent 
in SNe Ib/c and SNe lib systems. The difference is caused by 
the pre-SN orbital separations, which are, in general, smaller in 
SNe la compared to SN Ib/c systems. Simulations of SNe la in 
the literature have been calculated with small binary separations, 
assuming that the companion stars fill their Roche-lobe radii at 
the time of the explosions (Liu et al. 2012[|Pan et al.|[2012b|l. 


Nevertheless, one should keep in mind that the ejected stars that 
become HVSs originate from the closest pre-SN orbits ( |Tauris| 
2015|l, hence, SN-induced HVSs could suffer from SN ejecta 


impact effects. 


|Hirai & Yamada| ( |2015[ ) studied SN ejecta-companion inter¬ 
actions for iPTF 13bvn and found that the companion star may 
expand extremely after the SN impact. They considered binary 
systems with a companion mass of M 2 ~ 7.0 - 8.0 Mq and an 
orbital separation ratio of a IR 2 ~ 20 - 30. This result is dif¬ 
ferent from the predictions for these wide binary systems in our 
simulations. However, it is difficult to discuss the reasons that 
cause these differences because the descriptions of detailed hy- 
drodynamical impact simulations in |Hirai & Yamad^ ( |2015| l are 
limited. If the companion stars in binaries evolve to be red giant 
stars at the moment of CCSN explosion, the effect of the explo¬ 
sion is much greater because of the significantly reduced binding 


energy of the companion envelope ( |Hirai et al.|2014) l. However, 
as we have argued here, the probability that companion stars are 
red giant stars at the moment of a SN explosion is relatively low 
(< 15%). 

Finally, for hxed values of and Mg the calibrated fit¬ 
ting formulae for the removed mass or impact velocity in Ap¬ 
pendix A of |Tauris| ( |2015) l, for a 1.0 Mq MS star, are in good 
agreement with our hydrodynamical simulations for a 0.9 Mq 
MS star, yielding similar values for AM 2 and Vim. 


5. Summary and conclusions 

We have investigated the impact of SN ejecta on the compan¬ 
ion stars in CCSNe of Type Ib/c using the SPH code Stellar 
GADGET. We focused on 0.9 and 3.5 Mq MS companion stars 
constructed with the BEC stellar evolution code. Our main re¬ 
sults can be summarized as follows: 

1. The dependence of total removed mass (AM 2 ) and resulting 
impact velocity (vim) on the pre-SN binary separation (a) can 
be fitted with power-law functions, which is consistent with 
previous work on SNe la. The amount of accreted SN ejecta 
mass (AMacc) is found to be dependent on a with a power- 
law relationship. All three quantities are shown to decrease 
significantly with increasing a, as expected. 

2. The relation between AM 2 and the SN explosion energy (Eg) 
is linear. The quantities Vim and AMacc can better be htted 
with power-law functions and their values increase and de¬ 
crease, respectively, with larger values of Eej- 

3. With a hxed value of Eej, the amount of SN ejecta mass (Mej) 
is found to not affect AM 2 by much. However, more ejecta 
mass can be accumulated onto the companion star in simu¬ 
lations with larger value of Mej for a hxed value of Eej ■ The 
impact velocity, Vim, is sensitive to the momentum prohle of 
the outer SN ejecta. With our SN ejecta model, we hnd that 
the impact velocity, Vim, decreases at larger Mej. 

4. Our population synthesis calculations predict that in most 
CCSNe less than 5% of the MS companion mass can be re¬ 
moved by the SN impact (i.e., AM 2 /M 2 < 0.05). In addition, 
the companion star typically receives an impact velocity, Vim, 
of a few lOkms^*, and the amount of SN ejecta captured by 
the companion star after the explosion, AM^ec, is most often 
less than 10“^ Mq. 

5. Because a typical CCSN binary companion is relatively mas¬ 
sive and can be located at a large pre-SN distance, we do not 
expect, in general, that the effects of the SN explosion on the 
post-impact stellar evolution will be very dramatic. 

6. In the closest pre-SN systems, the MS companion stars 
are affected more strongly by the SN ejecta impact, lead¬ 
ing to AM 2 /M 2 - 0.10, Vim - 100 kms^* and AMacc - 
4 X 10“^ Mq, depending on the mass of the companion star. 
In addition, these stars are signihcantly bloated as a conse¬ 
quence of internal heating by the passing shock wave. 

7. It is possible that the SN-induced high velocity stars (HVSs), 
or more ordinary, less fast, runaway stars, may be sufficiently 
contaminated to be identihed by their chemical peculiarity as 
former companion stars to an exploding star if mixing pro¬ 
cesses are not efficient on a long timescale. 
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